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Abstract 

Both families of high T c superconductors, iron pnictides and cuprates, exhibit material dependence of superconductivity. 
Here, we study its origin within the spin fluctuation pairing theory based on multiorbital models that take into account 
realistic band structures. For pnictides, we show that the presence and absence of Fermi surface pockets is sensitive to 
the pnictogen height measured from the iron plane due to the multiorbital nature of the system, which is reflected to 
the nodeless/nodal form of the superconducting gap and T c . Surprisingly, even for the cuprates, which is conventionally 
modeled by a single orbital model, the multiorbital band structure is shown to play a crucial role in the material 
dependence of superconductivity. In fact, by adopting a two orbital model that considers the d z 2 orbital on top of 
the d x 2_ y 2 orbital, we can resolve a long standing puzzle of why the single layered Hg cuprate have much higher 
T c than the La cuprate. Interestingly, here again the apical oxygen height measured from the Cu02 plane plays an 
important role in determining the relative energy difference between d x i_ y i and d z 2 orbitals, thereby strongly affecting 
the superconductivity. 
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1. Introduction 

Superconductivity in the iron-based pnictide LaFeAsO 
doped with fluorine discovered by Hosono's groupfHQ and 
subsequent increase of the transition temperature (T c ) in 
the same family of compounds are remarkable as the first 
non-copper compound that has T c 's exceeding 50 K.Q 
Theoretical studies have shown that the electron-phonon 
Coupling is too weak to account for the high T C Q, and ex- 
ploring electronic pairing mechanisms has become a chal- 
lenge. 

An interesting observation of this series of material is 
their strong material dependence of superconductivity. T c 
ranges from 5K in LaFePOQ to 55K0] in SmFeAsO, 
where the importance of the lattice structure has been 
pointed out[6|. A number of experiments on arsenides 
suggest fully-open superconducting gap, consistent with 
the proposal of sign reversing s-wave pairing, 0, Q but 
experiments have shown presence of nodes in the super- 
conducting gap of LaFePOHEI .11], suggesting that even 
the gap is material dependent. 

The strong material dependence of superconductivity 
in the pnictides reminds us of the high T c cuprates, 





'Corresponding author. 
Email address: kurokiSvivace . e-one .uec . ac .up (Kazuhiko 
Kuroki) 



Fig.l Anion height measured from the cation plane in the 
iron pnictides (top) and the cuprates (bottom). 



which also exhibit strong material dependence of T c . It 
is well known that T c varies with the number of layers 
per unit cell, but an even more basic issue is the sig- 
nificant difference in T c within the single-layered mate- 
rials, i.e., La2- 2; (Sr/Ba) a ;Cu04 with a maximum T c of 
about 40K versus HgBa 2 Cu0 4+5 with a T c ~ 90K. Phe- 
nomenologically, it has been recognized that the mate- 
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rials with T c ~ 100K tend to have "round" Fermi sur- 
faces, while the Fermi surface of the "low T c " La system is 
closer to a square shape which implies a relatively better 
nesting [l~2l fl3|. 

Conventionally, the cuprates with a rounded Fermi sur- 
face have been modeled by a single-band model with large 
second (t 2 (> 0)) and third (ts(< 0)) neighbor hopping in- 
tegrals ((|t2| + I^D/I^il ~ 0.4), while the La system has 
been considered to have smaller t 2 ,t 3 ((\t 2 \ + |t3|)/|£i| ~ 
0.1). This, however, has brought about a contradiction 
between theories and experiments. Namely, while some 
phenomenological[3| and t-J model[l5|, [l6[ studies give 
a tendency consistent with the experiments, a number of 
many-body approaches for the Hubbard-type models with 
realistic values of on-site U show suppression of supercon- 



ductivity for large t 2 > and/or t 3 < 0jl7lll8l.ll 



In this paper, we show that the experimentally observed 
material/lattice structure dependence of superconductiv- 
ity in both pnictides [2oT| and cuprates 21 1 can be under- 
stood by analyzing the superconductivity within the spin 
fluctuation pairing theory that takes into account the real- 
istic band structure. A surprising finding is that not only 
the pnictides, but also the cuprates have band structures 
where multiple d orbitals play an important role. Inter- 
estingly, in both families, the anion height measured from 
the cation plane (Fig.l) turns out to be one of the key 
parameters that controls the multiorbital band structure 
and thus the superconductivity. These studies show that a 
combination of effective model construction based on first 
principles band calculation and the application of many 
body theory can give a realistic description on the material 
dependence of superconductivity in correlated systems. 



2. Construction of multiorbital models 

2.1. Pnictides 

We start with the band structure of LaFeAsO. LaFeAsO 
has a layered structure, where Fe atoms form a square lat- 
tice in each layer, which is sandwiched by As atoms. Due 
to the tetrahedral coordination of As atoms, there are two 
Fe atoms per unit cell. The experimentally determined 
lattice constants are a = 4.03A and c = 8.74A, with two 
internal coordinates ZLa = 0.142 and zas = 0.651. [22j We 
have first obtained the band structure with plane-wave 
basis 2J| , which is then used to construct the maximally 
localized Wannier functions 23] . These Wannier functions 
have five orbital symmetries {dj,z 2 -R 2 , dxz, dyz, dx 2 -Y 2 , 
dxY j where X, Y, Z refer to those for this unit cell with two 
Fe sites as shown in Fig. 2(a). The two Wannier orbitals 
in each unit cell are equivalent in that each Fe atom has 
the same local arrangement of other atoms. We can then 
take a unit cell that contains only one orbital per symme- 
try by unfolding the Brillouin zone, and we end up with 
an effective five-band model on a square lattice, where x 
and y axes are rotated by 45 degrees from X-Y. We refer 
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Fig. 2 (a) The original (dashed lines) and reduced (solid) 
unit cells with • (Fe), V (As below the Fe plane) and 
A (above Fe). (b) The band structure (left) of the five- 
band model for LaFeAsO, and the Fermi surface (right) at 
k z = for n = 6.1. The main orbital characters of some 
portions of the bands and the Fermi surface are indicated. 
The dashed horizontal line in the band structure indicates 
the Fermi level for n — 6.1. The short arrow in the band 
structure indicates the position of the Dirac cone closest 
to the Fermi level. The gray areas in the Fermi surface 
around the zone corners represent the 7 Fermi surface, (c) 
The portion of the band that has mainly the g?x 2 -y 2 (solid 
red) and d z i (dashed blue) orbital character. 



all the wave vectors in the unfolded Brillouin zone here- 
after. We define the band filling n as the number of elec- 
trons/number of sites (e.g., n = 10 for full filling). The 
doping level x in LaFeAsOi-^Fj, is related to the band 
filling as n = 6 + x. 

The five bands arc heavily entangled as shown in 
Fig. 2(b) (left) reflecting strong hybridization of the five 
3d orbitals, which suggests that the minimal model for the 
pnictides should include all the five orbitals. In Fig. 2(b) 
(right), the Fermi surface for n = 6.1 (corresponding to 
x = 0.1) obtained by ignoring the inter- layer hoppings 
is shown in the two-dimensional unfolded Brillouin zone. 
The Fermi surface consists of four pieces: two concentric 
hole pockets (denoted here as ai, a 2 ) centered around 
(k x , k y ) — (0, 0), two electron pockets around (tt, 0) (f5\) or 
(0,7r) (/3 2 ), respectively, a, (/%) corresponds to the Fermi 
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surface around the TZ line (MA in the original Brillouin 
zone) in the first-principles band calculation. [25| Besides 
these pieces of the Fermi surface, there is a portion of the 
band near (tt, tt) that touches the Ep, so that the portion 
acts as a "quasi Fermi surface (7)" around (tt, tt), which 
has in fact an important contribution to the spin suscep- 
tibility as we shall see below. As for the orbital character, 
a and portions of /3 near Brillouin zone edge have mainly 
dxz and dyz character, while the portions of j3 away from 
the Brillouin zone edge and 7 have mainly dx 2 -Y 2 orbital 
character (see also Fig. 2(c)). 

To be more precise, the above band structure is sen- 
sitive to the lattice structure. In the upper two panels 
of Fig. 3, we compare the band structure of LaFePO and 
NdFcAsO. A large difference between the two materials 
lies in the band structure near the wave vector (tt,tt). In 
LaFePO, the 7 Fermi surface originating from the d X 2 -Y 2 
band around (tt, tt) largely sinks below the Fermi level, and 
in turn the dz 2 band, which is not effective for supercon- 
ductivity, rises up. On the other hand, in NdFeAsO, the 
dx 2 -Y 2 band rises up even more than in LaFeAsO. The 
origin of this band structure variation is mainly due to 
the variation of the pnictogen height hp n 2(], 2a, j 
which varies from hp n = 1.14A in LaFePO to hpj, 



mm, 

1.38A 

in NdFeAsO. [6| In fact, in the lower two panels of Fig. 3, 
we show band structures of LaFeAsO with hypothetical 
lattice structures, where we fix the lattice constants and 
change only the height to those of NdFeAsO or LaFePO. 
We can see that the band structure around (tt, tt), namely, 
the position of the dx 2 -Y 2 and dz 2 bands are determined 
by the pnictogen height. 

The origin of this height sensitivity of the band struc- 
ture is the following. In Fig. 2(c), we plot the d X 2 ~y 2 an d 
d Z 2 portions of the bands. Around (tt, tt), the lower por- 
tion of dx 2 -Y 2 and the upper portion of dz 2 lies close 
to the Ep. When the pnictogen height becomes low, the 
Wannier orbitals widely spread toward the pnictogen, so 
that the band width tends to become large. Therefore, the 
dx 2 -Y 2 portion around (tt,tt) sinks below Ep, while the 
dz 2 portion rises up. 
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Fig. 3 From top to bottom :The band structure and the 
Fermi surfaces of NdFeAsO, LaFePO, hypothetical struc- 
ture of LaFeAsO with h Pn = 1.38A and 1.14A. 



2.2. Cuprates 

In contrast to the pnictides, theoretical studies on 
cuprates has mostly been done for single band models, 
or the three band model that takes into account the oxy- 
gen orbitals that strongly hybridizes with the Cu d x 2_ y 2 
orbital. As mentioned in the introduction, within these 
single orbital approaches, it is difficult to understand the 
experimental observation that the more Fermi surface is 
rounded, higher the T c . To resolve the discrepancy be- 
tween experimental observations and theory, we consider 
a two- orbital model, where we take into account the d z 2 
orbital on top of the d x 2_ y 2 21|. In fact, for the La sys- 



tem, it has long been known that a band with a stron g d 7 2 



character lies rather close to the Fermi energy 28, 29, 30]. 



More recently, it has been discussed in refs. 31 , 12j that 



the shape of the Fermi surface is determined by the en- 
ergy level of the "axial state" consisting of a mixture of 
Cu d z 2-0 p z and Cu 4s orbitals, and that the strength 
of the d z 2 contribution causes the difference in the Fermi 
surface shape between the La and Hg systems. Namely, 
the d z 2 contribution is large in the La system making the 
Fermi surface closer to a square, while the contribution 
is small in the Hg system making the Fermi surface more 
rounded. In Fig. 4, we show the first-principles [24| result 
for band structures in the two-orbital model for the La and 
Hg systems, obtained by constructing maximally local- 
ized Wannier orbitals (23|. The lattice parameters adopted 
here are experimentally determined ones for the doped 
materials (32T l33j . We can confirm that in the La system 
the main band (usually considered as the ll d x 2_ y 2 band") 
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has in fact a strong d z 2 character on the Fermi surface 
near the N point, which corresponds to the wave vectors 
(tt, 0), (0, 7r) in the Brillouin zone of the square lattice. The 
d z 2 contribution is seen to "push up" the van Hove singu- 
larity (vHS) of the main band, resulting in a seemingly 
well nested (square shaped) Fermi surface. In the Hg sys- 
tem, on the other hand, the d z 2 band stays well away from 
Ep, and consequently the vHS is lowered, resulting in a 
rounded Fermi surface. 

If we estimate in the two-orbital model the ratio ( | i 2 1 + 
|i3|)/|ii| within the d x 2_ y 2 orbitals, we get 0.35 for the 
La system against 0.41 for Hg, which are rather close to 
each other. This contrasts with the situation in which the 
model is constrained into a single band model that con- 
siders one kind of Wannier orbital to reproduce the main 
band that intersects the Fermi energy. There, the Wannier 
orbital has mainly d x 2_ y 2 character, but has "tails" with 
a d z 2 character especially for the La system. Then the ra- 
tio (1^1 + |*3|)/|ii| in the single-orbital model reduces to 
0.14 for La against 0.37 for Hg, which is just the conven- 
tional view mentioned in the Introduction. From this, we 
can confirm that it is the d z 2 contribution that makes the 
Fermi surface in the La system square shaped, while the 
"intrinsic" Fermi surface of the high T c cuprate family is, 
as in the Hg system, rounded (which is actually due to the 
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Fig. 4 The band structure of the two (d x 2_ y 2-d z 2) orbital 
model for La 2 Cu0 4 (left) and HgBa 2 Cu0 4 (right). The 
top (middle) panels depict the strength of the d x 2_ y 2 (d x n) 
characters with thickened lines, while the bottom panels 
the Fermi surfaces (for a total band filling n = 2.85). 



Cu 4s orbital as shown in Ref. 3l|, [12J , whose contribution 
is effectively considered in the present model). 

3. Many-body Hamiltonian 

For the many body part of the Hamiltonian, we con- 
sider the standard interaction terms that comprise the 
intra-orbital Coulomb U, the inter-orbital Coulomb U', 
the Hund's coupling J and the pair-hopping J'. The many 
body Hamiltonian then reads 



H 



IJ ftU a 



i 

E 



,(1) 



where i, j denote the sites and /i, v the orbitals, and tfj is 
the transfer energy obtained from the maximally localized 
Wannier orbitals. 

For these models, we apply spin fluctuation theory and 
solve the Eliashberg equation to analyze superconductiv- 
ity. Multiorbital random phase approximation (RPA) is 
described in e.g. Refs.[3~4ll35|. In the present case, Green's 
function Gi m (k) (k = (k,ioj n j) is a m x rn matrix, where 
m is the number of orbitals. The irreducible susceptibility 
matrix 

xl,i M (q) = -J2 G ^( k + l) G Ui2( k ) (2) 

k 

(h = 1, m) has m 2 x m 2 components, and the spin and 
the charge (orbital) susceptibility matrices are obtained 
from matrix equations, 



Xs(q) = 
Xc(q) = 



x°(q) 



l-S X °(q) 

X°(q) 
l + C X °(q) 



where 




(3) 
(4) 



h = h = h = h 

h = h h = h 

h = h 7^ h = h 

h — U 7^ h = h 



The Green's function and the effective singlet pairing in- 
teraction, 

V s (q) = ls X s (q)S \c X c (q)C + \(S + C), (5) 
are plugged into the linearized Eliashberg equation, 

q 22*3*5 £'6 

Gi 2 i 5 { k ~ q)4>hi 6 {k - q)Gi 3 i 6 (q - fc)(6) 
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Here, T c is the temperature where the eigenvalue A reaches 
unity, so when A is calculated at a fixed temperature for 
different situations/materials, the value can be considered 
as a qualitative measure for the T c . 



4. Anion height dependence of superconductivity 

4--1- Pnictides 

We first look into the superconductivity of pnictides. We 
apply RPA to the 5 orbital model, where we take U = 1.2, 
U' = 0.9, J = J' = 0.15eV, the temperature is fixed at 
T = 0.02 eV and 32 x 32 x 4 fc-point meshes and 512 
Matsubara frequencies are taken. Let us first look at the 
the spin susceptibility of LaFcAsO. In the top panels of 
Fig. 5, we show Xs3333 and Xs4444j which represent the spin 
correlation within dyz and dx 2 -Y 2 orbitals, respectively. 
Xs3333 peaks solely around (jr, 0) and (0,7r), which reflects 
the nesting between the XZ, FZ-charactered portions of a 
and P Fermi pockets as shown in a bottom panel of Fig. 5. 
On the other hand, x$4444 has peaks around (tt, 0), (0,7r) 
and around (tt,tt/2)/(tt/2,tt) as well. The former is due 
to the nesting between the 7 quasi Fermi surface and the 
d x z-Y 2 portion of the (3 Fermi surface as was first pointed 
out in ref. 7], while the latter originates from the nest- 
ing between the dx 2 -Y 2 portion of the fi\ and /3a Fermi 
surfaces. Q The (tt, 0), (0, tt) feature is consistent with the 
stripe (i.e., collinear) antiferromagnetic order for the un- 
doped case,which was suggested by transport and optical 
reflectance [361], and further confirmed by neutron scatter- 
ing experiments. [22| The stabilization of such an antiferro- 
magnetic ordering has also been pointed out in first prin- 
ciples calculations. 37, 36, |38| Furthermore, the presence 
of spin fluctuations near the wavevector (tt, 0), (0, tt) in the 
unfolded Brillouin zone has indeed been confirmed in an 
inelastic neutron scattering experiment [39j| . 

When the pairing is mediated by spin fluctuations, the 
gap on the Fermi surface tends to change its sign across the 
wave vector at which the spin fluctuations develop. Since 
there are multiple spin fluctuation modes in our model, 
the superconducting gap should be determined by the co- 
operation or competition between these multiple modes. 
Specifically, the a-f3 and 7-/3 nestings tend to favor the 
fully-gapped, sign-reversing s-wave (s±-wave), in which 
the gap changes sign between a and /3 and also between 
7 and /3, but has a constant sign on each pocket as shown 
in Fig.5.0. This kind of gap has in fact been considered 
in the context of raising T c in spin fluctuation mediated 
pairing in that the nodes of the gap, which is a necessary 
evil for spin fluctuation mediat ed p airing, do not have to 
intersect the Fermi surface 140 



41] . On the other hand, 




/3i-/?2 nesting tends to change the sign of the gap between 
these pockets, which can result in either d-wave or an s- 
wave pairing with nodes on the /3 Fermi surface, as shown 
in Fig.5 HHU. 

When the pnictogen atom is at low positions, the /3 — 7 
nesting is not effective, and the competition between 



Fig.5 Top panel: RPA result of the diagonal compo- 
nent of the spin susceptibility matrix, x S 3333(left) and 
X S 4444(right) (3 : YZ, 4 : X 2 - Y 2 ) for LaFeAsO. Bottom 
panel: Nesting is shown for the Fermi surface for orbitals 
XZ, YZ (left) and for X 2 - Y 2 (right). Here the thickness 
of the Fermi surface represents the strength of the respec- 
tive orbital character. The gray areas around the corners 
in the left panel indicates the 7 "quasi" Fermi surface. 



(7r/0)/(0,7r) and (tt, 7r/2)/(7r/2, tt) spin fluctuations results 
in a frustration in momentum space, where the gap tends 
to change the sign of the gap across these vectors. Conse- 
quently, the s-wave gap has nodes intersecting the /3 Fermi 
surface, and this pairing is nearly degenerate with cZ-wave 
pairing. When the pnictogen is high, the /3-7 nesting be- 
comes more effective than the /?i-/?2 nesting, and the two 
kinds of of (tt/0)/ (0, tt) spin fluctuations, those originating 
from dxz/dyz and dx 2 -Y 2 orbitals, cooperates without 
experiencing frustration to result in the fully gapped si- 
wave pairing. 

The absence/presence of the frustration affects also the 
strength of the superconducting instability, where the fully 
gapped case has a stronger tendency toward pairing. This 
can be seen from the eigenvalue of the Eliashberg equation 
plotted as a function of the pnictogen height in the bot- 
tom panel of Fig. 6, where we adopt the 5 band model for 
LaFeAsO, but vary the height. Thus the pnictogen height 
hp n acts as a switch between higher T c nodeless and low 
T c nodal pairings, which explains the difference between 
NdFeAsO and LaFePO in both T c and the gap. 

Here, we have focused on the effect of the 7 Fermi sur- 
face around (tt, tt) in the context of the lattice structure 
dependence, but its effect has recently been studied in de- 
tail from the viewpoint of the electron-electron interac- 
tion and the band filling 0, |44|. Also, a real space pic- 
ture for the presence/absence of the nodes in the gap has 
been discussed 45] . The electron correlation effects beyond 
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Fig. 6 Top and middle panels: the fully-gapped s± wave, 
the nodal s± wave and d wave gap are schematically 
shown. The solid red (dashed blue) curves represent pos- 
itive (negative) sign of the gap. The arrows indicate the 
dominating nesting vectors. Bottom: the s-wave and d- 
wave eigenvalues of the Eliashberg equation plotted as 
functions of As height for the hypothetical lattice struc- 
ture of LaFeAsO. The figures in the left and right insets 
show schematically the presence or absence of frustration. 



RPA has also been taken into account in the fluctuation 
exchange (FLEX) study 05^ and a functional renormal- 
ization group study [13], which also finds a similar effect of 
the 7 Fermi surface on the presence/absence of nodes in 
the superconducting gap. 

The above scenario is consistent with the presence of 
line nodes in the superconducting gap of the low T c phos- 
phides 0, [HI U | ■ On the other hand, the fully gapped 
s±-wave is difficult to detect experimentally because the 
absolute value of the gap is essentially the same as that 
of a simple s-wave gap, with no characteristic angular de- 
pendence. Nevertheless, a recent phase sensitive STM ex- 
periment has clearly detected the sign reverse between the 
Fermi surfaces separated by (ir, 0), (0, 7r) [48] . 



4-2. Cuprates 

We now move on to the superconductivity in the 
cuprates. For the electron-electron interactions, we take 
U = 3.0eV, J — J' — 0.3eV, which gives the interor- 
bital U' = U - 2J = 2.4eV. The temperature is fixed 



at ksT = O.OleV. As for the band filling (number of elec- 
trons/site), we concentrate on the total n = 2.85, for which 
the main band has 0.85. Here we apply the multiorbital 
FLEX[34ll!{|[ioj]> which takes into account the self-energy 
correction to the Green's function self-consistently, for the 
three-dimensional lattice taking 32 x 32 x 4 fc-point meshes 
and 1024 Matsubara frequencies. 

Here again, we focus on how the anion height, namely, 
the apical oxygen height measured from the Cu02 plane 
affects the band structure and hence superconductivity. 
This is motivated by the fact that the energy level offset 
between d x 2_ y 2 and d z 2 orbitals should be controlled (at 
least partially) by the ligand field, hence by the height, 
ho, of the apical oxygen [2 9]. To single out this effect, let 
us examine the two-orbital model of the La cuprate for 
which we increase ho from its original value 2.4lA with 
other lattice parameters fixed. In Fig. 7(a), we plot the 
Eliashberg equation eigenvalue of the d-wave superconduc- 
tivity as a function of ho- We can see that A monoton- 
ically increases with the height. As seen from the inset 
of Fig. 7(b), AE, defined as the difference of the on-site 
energy between d x i_ y i and d z 2 orbitals, is positively cor- 
related with ho as expected, and importantly, Fig. 7(b) 
shows that the increase in A is positively correlated with 
the increase in AE. Hence the superconductivity turns out 
to be enhanced as the d z 2 band moves away from the main 
band. Note that this occurs despite the Fermi surface be- 
coming more rounded with larger AE, namely, the effect of 
the orbital character (smaller d z 2 contribution) dominates 
over the Fermi surface shape effect. Conversely, the strong 
d z 2 orbital character in the Fermi surface around the wave 
vectors (tt, 0), (0, tt) works destructively against d-wave su- 
perconductivity. Physically, this may be understood as 
follows, d-wave superconductivity occurs due to pair scat- 
tering from ~ (tt, 0) to ~ (0,7r) and vice versa mediated 
by antiferromagnetic spin fluctuations (Fig. 7(c)). When 
the states around (w, 0), (0, tt) has strong d z 2 character, 
superconductivity is suppressed since d-wave pairing has a 
rough tendency for higher T c in bands that are nearly half 
filled, whereas the d z 2 orbital here is nearly full filled (has 
holes only around (tt, 0), (0,7r)). 

In Fig. 7, we have also plotted the corresponding values 
of the Hg system obtained with the actual lattice struc- 
ture. We can see that AE in Hg is indeed larger than 
that in La as expected, but actually AE ~ 2.2eV for Hg 
is larger than AE ~ 1.3eV, which is the value the La sys- 
tem would take for ho — 2.8A. Consequently, A for Hg 
is somewhat larger than that for the La system with the 
same value of ho- This implies that there are also some 
effects other than the apical oxygen height that also en- 
hance AE (i.e., lower the d z 2 level with respect to the 
d x 2_ y 2 level) in the Hg system, thereby further favoring 
d-wave superconductivity. In this context, the present re- 
sult reminds us of the so-called "Maekawa's plot" , where 
a positive correlation between T c and the level of the api- 
cal oxygen p z hole was pointed out [51] . Since a higher p z 
hole level (i.e., a lower p z electron level) is likely to lower 
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Fig. 7 The eigenvalue of the Eliashberg equation A (red cir- 
cles) when ho(&) or AE(b) is varied hypothetically in the 
lattice structure of La2Cu04. Blue diamond indicates the 
eigenvalue of HgBa2Cu04. Inset in (b) shows the relation 
between ho and AE. In (c), the pair scattering of d z 2 
holes around (tt, 0) (0, tt) is illustrated. 



interesting to note that in both families, the anion height 
measured from the cation plane plays an important role, 
although its effect on the band structure appears in a dif- 
ferent manner : in the cuprates the the energy difference 
between d x 2_ y 2 and d z 2 is affected by the ligand field, while 
in the pnictides, the relative position between the dx 2 -Y 2 
and dz 2 bands are affected mainly due to the variation of 
the band width, namely, the former (latter) band lies in 
the upper (lower) portion of the band structure, so that 
when the band width becomes wider upon lowering the 
height, the location of the lower band edge of dx 2 -Y 2 and 
that of the upper band edge of dz 2 (both around the wave 
vector (it, tt) )is exchanged. The good description of the 
material dependence of the superconductivity within the 
present approach suggests commonality between the two 
families in that the spin degrees of freedom are playin g an 
important role in the occurrence of superconductivity [52|, 
although there are differences in the strength of the elec- 
tron correlation. 

The present studies show that a combination of effec- 
tive model construction based on first principles band 
calculation and application of many body theory can 
give a realistic description on the material dependence 
of superconductivity in correlated systems, and may 
even suggest a possibility of theoretically predicting new 
superconductors in the near future. 
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E z 2 because the d z 2 Wannier orbital in our model consists 
mainly of Cu3c^2 and apical oxygen 2p z orbitals, the pos- 
itive correlation between AE and T c found here is indeed 
consistent with Maekawa's plot. 

5. Conclusion 

In this paper, we have focused on the origin of the mate- 
rial dependence of superconductivity in the two families of 
high T c superconductors, pnictides and cuprates. We have 
analyzed the superconductivity within the spin fluctuation 
pairing theory applied to models that take into account 
the realistic multiorbital band structure. It is shown that 
not only in the iron pnictides, but even in the cuprates, 
the multiorbital band structure is affected by the lattice 
structure, which in turn affects superconductivity. It is 
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